Applications of random matrix theory 
to condensed matter and optical physics 



C.W.J. Beenakker 

Instituut-Lorentz, Universiteit Leiden, P.O. Box 9506, 2300 RA Leiden, The Netherlands 

Two chapters for The Oxford Handbook of Random Matrix Theory, edited by G. Akemann, 
J. Baik, and P. Di Francesco, to be published by Oxford University Press. 



Contents 



I Condensed Matter Physics 



I. Introduction Q] 

II. Quantum wires [3 

A. Conductance fluctuations 

B. Nonlogarithmic eigenvalue repulsion 

C. Sub-Poissonian shot noise [5| 

III. Quantum dots [B| 

A. Level and wave function statistics 

B. Scattering matrix ensembles [B| 

C. Conductance distribution @ 

D. Sub-Poissonian shot noise 

E. Thermopower distribution 

F. Quantum-to-classical transition @ 

IV. Superconductors 

A. Proximity effect 

B. Gap fluctuations 

C. From mesoscopic to microscopic gap 

D. Quantum-to-classical transition 



II Classical and Quantum Optics 

V. Introduction [ll] 

VI. Classical optics [ll] 

A. Optical speckle and coherent backscattering [ll] 

B. Reflection from an absorbing random medium [l^ 

C. Long-range wave function correlations [lj| 

D. Open transmission channels [l4| 



VII. Quantum optics [lCj 

A. Grey-body radiation [if] 

B. RMT of a chaotic laser cavity [TtI 

Acknowledgments [Igl 

References 03 



Part I 

Condensed Matter 
Physics 

I. INTRODUCTION 

Applications of random matrix theory (RMT) to con- 
densed matter physics search for universal features in the 
electronic properties of metals, originating from the uni- 
versality of eigenvalue repulsion. Eigenvalue repulsion is 
universal, because the Jacobian 



j = n \Ej - Ei 



(1.1) 



%<3 



of the transformation from matrix space to eigenvalue 
space depends on the symmetry of the random matrix 
ensemble (expressed by the index (3 £ {1,2,4}) — but is 
independent of micros copic pr operties such as the mean 
eigenvalue separation ( Meh67l ). This universality is at 
the orig in of the rema rkable success of RMT in nuclear 
physics (|Bro8lUWei09h . 

In condensed matter physics, the applications of RMT 
fall into two broad categories. In the first category, one 
studies thermodynamic properties of closed systems, such 
as metal grains or semiconductor quantum dots. The 
random matrix is the Hamiltonian H. In the second cat- 
egory, one studies transport properties of open systems, 
such as metal wires or quantum dots with point con- 
tacts. Now the random matrix is the scattering matrix 
S (or a submatrix, the transmission matrix t). Appli- 
cations in both categories have flourished with the de- 
velopment of nanotechnology. Confinement of electrons 
on the nanoscale in wire geometries (quantum wires) and 
box geometries (quantum dots) preserves their phase co- 
herence, which is needed for RMT to be applicable. 

The range of electronic properties addressed by RMT 
is quite broad. The selection of topics presented in this 
Chapter is guided by the desire to show those applica- 
tions of RMT that have actually made an impact on ex- 
periments. For a more complete coverage of topics and 
a more compre hensive list of reference s we suggest a few 
review articles (lBee97l : lGuh98l lAlhOd h 
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FIG. 1 Fluctuations as a function of perpendicular magnetic 
field of the conductance of a 310 nm long and 25 nm wide Au 
wire at 10 mK. The trace appears random, but is completely 
reproducible from one measurement to the next. The root- 
mean-square of the fluctuations is 0.3 e 2 /h, which is not far 
from the theoretical result ^1/15 e 2 /h [Eq. (JUJ with (3 = 2 
due to the magnetic field and a reduced conductance quantum 
Go = e 2 /h due to the str ong spin-orbit scattering in Au]. 
Adapted from Ref. (|Was86t) , 

II. QUANTUM WIRES 

A. Conductance fluctuations 

In the 1960's, Wigner, Dyson, Mehta, and others dis- 
covered that the fluctuations in the energy level density 
are governe d by lev el repulsion and therefore take a uni- 
versal form (Por65). The universality of the leve l fluctu- 
ations is expressed by the Dyson-Mehta formula (|Dvs63h 
for the variance of a linear statistic 1 A = J2 n a (E n ) on 
the energy levels E n . The Dyson-Mehta formula reads 



Var A 



1 1 

7? 



dk\a(k)\ 2 k, 



(2.1) 



where a(k) — dE e lkE a(E) is the Fourier transform 
of a(E). Eq. (|2.1| shows that: 1. The variance is inde- 
pendent of microscopic parameters; 2. The variance has 
a universal 1 //3-dependence on the symmetry index. 

In a pair of seminal 1986-papers (llmr86l : lAlt86h . Imry 
and Altshuler and Shklovkskh proposed to apply RMT 
to the phenomenon of universal conductance fluctuations 
(UCF) in metals, which was discovered using diagram- 
matic pertu rbation theory by Altshuler ( Alt85t ) and Lee 
and Stone (|Lee85t i. UCF is the occurrence of sample- 
to-samplc fluctuations in the conductance which are of 
order e 2 /h at zero temperature, independent of the size 



The quantity A is called a linear statistic because products of 
different E n 's do not appear, but the function a(E) may well 
depend non-linearly on E. 



of the sample or the degree of disorder — as long as the 
conductor remains in the diffusive metallic regime (size L 
large compared to the mean free path I, but small com- 
pared to the localization length £) . An example is shown 
in Fig. [Q 

The similarity between the statistics of energy levels 
measured in nuclear reactions on the one hand, and the 
statistics of conductance fluctuations measured in trans- 
port exp_eriinentson the other hand, was used by Stone et 
al. (|Mut87l : ISto9l"ri to construct a random matrix theory 
of quantum transport in metal wires. The random ma- 
trix is now not the Hamiltonian H , but the transmission 
matrix t, which determines the conductance through the 
Landauer formula 



G=G Trttt = G ^T„. 



(2.2) 



The conductance quantum is Go = 2e 2 //i, with a fac- 
tor of two to account for spin degeneracy. Instead of 
repulsion of energy levels, one now has repulsion of the 
transmission eigenvalues T n , which are the eigenvalues of 
the transmission matrix product ttf . In a wire of cross- 
sectional area A and Fermi wave length Xp, there are of 
order N ~ A/Xp propagating modes, so t has dimen- 
sion N x N and there are N transmission eigenvalues. 
The phenomenon of UCF applies to the regime N 3> 1, 
typical for metal wires. 

Random matrix theory is based on the fundamen- 
tal assumption that all correlations between the eigen- 
values are due to the Jacobian J = J\i<j l^i — 
from matrix elements to eigenvalues. If all correlations 
are due to the Jacobian, then the probability distribu- 
tion P(T 1 ,T 2 , . . . T N ) of the T„'s should have the form 
P cx JY\iP{Ti), or equivalently, 

P({T„})cx ex V \-(3{y j u{T i ,T j )+Y J V{T i ))\, (2.3) 



i<3 



(2.4) 



with V = -/3' 1 Inp. Eq. ([23)1 has the form of a Gibbs 
distribution at temperature for a fictitious system 
of classical particles on a line in an external potential V, 
with a logarithmically repulsive interaction u. All mi- 
croscopic parameters are contained in the single function 
V(T). The logarithmic repulsion is independent of mi- 
croscopic parameters, because of its geometric origin. 

Unlike the RMT of energy levels, the correlation func- 
tion of the TVs is not translationally invariant, due to 
the constraint < T n < 1 imposed by unitarity of the 
scattering matrix. Because of this constraint, the Dyson- 
Meht a formu la (|2.ip needs to be modified, as shown in 
Ref. (|Bee93aT) . In the large-A^ limit, the variance of a 
linear statistic A = J2 n f{T n ) on the transmission eigen- 
values is given by 



1 1 f°° 

VarA=-— dfc|F(fc)| 2 fctanh(7rfc). 
P J 



(2.5) 
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The function F(k) is defined in terms of the function 
f(T) by the transform 



F(k) 



dxe ikx f 



(2.6) 



The formula (|2.5p demonstrates that the universality 
which was the hallmark of UCF is generic for a whole 
class of transport properties, viz. those which are lin- 
ear statistics on t he tran smission eigenvalues. Examples, 
reviewed in Ref. (|Bee97h . are the critical-current fluctu- 
ations in Josephson junctions, conductance fluctuations 
at normal-superconductor interfaces, and fluctuations in 
the shot-noise power of metals. 



B. Nonlogarithmic eigenvalue repulsion 



The probability distribution (|2.3|) was justified by a 
maximum-e ntropy principle for an ensemble of quasi-lD 
conductors (|Mut87l ; ISto9lh . Quasi-lD refers to a wire 
geometry (length L much greater than width W). In 
such a geometry one can assume that the distribution of 
scattering matrices in an ensemble with different realiza- 
tions of the disorder is only a function of the transmission 
eigenvalues (isotropy assumption) . The distribution (|2.3[) 
then maximizes the information entropy subject to the 
constraint of a given density of eigenvalues. The function 
V(T) is determined by this constraint and is not specified 
by RMT. 

It was initially believed that Eq. (|2.3[) would provide an 
exact description in the qua si- ID l imit L 3> W, if only 
V(T) were suitably chosen l Sto9ll ). However, the gen- 
eralized Dyson-Mehta formula (|2.5[) d emonstra tes that 
RMT is not exact in a quantum wire (|Bee93ah . If one 
computes from Eq. (|2.5p the variance of the conductance 
QUI) [by substituting f[T) = G T], one finds 



the T's? Or is the true interaction a many-body inter- 
action, which cannot be reduced to the sum of pairwise 
interactions? This transport problem has a counterpart 
in a closed system. The RMT of the statistics of the 
eigenvalues of a random Hamiltonian yields a probability 
distribution of the form (|2.3p . with a logarithmic repul- 
sion bet ween t he energy levels (|Meh67h . It was shown b 
Efetov (|Efe83l ) and by Altshuler and Shklovskii (|Alt8' 
that the logarithmic repulsion in a disordered metal grain 
holds for energy separations small compared to the in- 
verse ergodic time h/T erg . 2 For larger s eparat ions the 
interaction potential decays algebraically ( Jal93l ). 

The way in which the RM T of qua ntum transport 
breaks down is quite different ( Bee93rJ ). The probabil- 
ity distribution of the transmission eigenvalues does in- 
deed take the form (|2.3p of a Gibbs distribution with 
a parameter-independent two-body interaction u(Ti,Tj), 
as predicted by RMT. However, the interaction differs 
from the logarithmic repulsion (|2.4p of RMT. Instead, it 
is given by 



u(Ti,Tj 



1 1*12} 



T,\ 



with Tn 



1 / cosh 2 : 



(2.9) 



The eigenvalue interaction (|2.9p is different for weakly 
and for strongly transmitting scattering channels: u — > 
— ln|Tj — Tj| for Ti,Tj — > 1, but u — > — -| In |Tj- — T,-| 
for Ti,Tj <C 1. For weakly transmitting channels it is 
twice as small as predicted by considerations based solely 
on the Jacobian, which turn out to apply only to the 
strongly transmitting channels. 

The nonlogarithmic interaction modifies the Dyson- 
Mehta formula for the variance of a linear stat istic. In- 
stead of Eq. (|2.5p one now has (|Bee93rJ : ICha93T ) 



VarG/Go = -0 



(2.7) 



independent of the form of V(T) . The diagrammatic 
perturbation theory ( Alt85l : ILee85l ) of UCF gives instead 



1 1 f°° 
VarA = - — / dk 



k\F{k)\< 



P 71-2 Jo 1 + cotanh(i7rfc) ' 



F(k) = f 

J — i 



dxe lkx f 



cosh 2 x 



(2.10) 
(2.11) 



VarG/Go = — /T 1 
15 



(2.8) 



for a quasi-lD conductor. The difference between the 
coefficients | and ^ is tiny, but it has the fundamen- 
tal implication that the interaction between the T's is 
not precisely logarithmic, or in other words, that there 
exist correlations between the transmission ei genvalue s 
over and above those induced by the Jacobian ( Bee93al ). 



The | — discrepancy raised the question what the 
true eigenvalue interaction would be in quasi-lD conduc- 
tors. Is there perhaps a cutoff for large separation of 



Substitution of f(T) = T now yields ^ instead of | for 
the coefficient of the UCF, thus resolving the discrepancy 
between Eqs. dH7J) and ([2T5]) . 

The result (|2.9p follows from the solution of a differen- 
tial equation which determines how the probability distri- 
bution of the T n 's changes when the length L of the wire 
is incremented. This differential equation has the form 
of a multivariate drift-diffusion equation (with L playing 
the role of time) for N classical particles at coordinates 



x n = arcoshT r 



-1/2 



The drift-diffusion equation, 



4 



2 The ergodic time is the time needed for a particle to explore the 
available phase space in a closed system. In a disordered metal 
grain of size L and diffusion constant D, one has r erg ~ L 2 /D. 



If the motion is ballistic (with velocity vp) rather than diffusive, 
one has instead r crg ~ L/vp. 



n— 1 

N N N 

= — In | sinh 2 x j — sinh 2 Xi | — — In | sinh 2a;, | , 



i=l j=i+l 



(2.12) 
(2.13) 



is known as the DMPK equation, after t he scientists wh o 
first studied its properties in the 1980's (|Dor82l ; [Mel88cJ) . 
(The equation itself appeared a decade earlier ( Bur72l).') 
The D MPK equation can be solved exactly (|Bee93bl : 
ICas95t) . providing the nonlogarithmic repulsion (|2.9p . 



C. Sub-Poissonian shot noise 

The average transmission probability T — 1/ L for dif- 
fusion through a wire is the ratio of mean free path I and 
wire length L. This average is not representative for a 
single transmission eigenvalue, because eigenvalue repul- 
sion prevents the TVs from having a narrow distribution 
around f. The eigenvalue density p(T) = (E n S(T -T n )) 
can be cal culated from th e DMPK equation (|2.12|) , with 
the result (|Dor84 lMel89t ) 



P{T) 



Nl 



1 



2L T\f\~^T 



for T min < T < 1, 



(2.14) 



in the diffusive metallic regime 3 I <C L <C £. The 
lower limit T m i n is determined by the normalization, 
Jp 1 dT p{T) = N, giving T min w 4e _i / 2/ with exponen- 
tial accuracy. 

The transmission eigenvalue density is bimodal, with a 
peak at unit transmission (open channels) and a peak at 
exponentially small transmission (closed channels). This 
bimodal distribution cannot be observed in the conduc- 
tance G cx ^2 n T n , which would be the same if all T n 's 
would cluster near the average T. The shot noise power 
iS oc En ^n(l — T n ) (the second moment of the time de- 
pendent current fluctuations) provides more information. 

The ratio of shot noise power and conductance, defined 
in dimensionless form by the Fano factor 



F = 



(2.15) 



3 The localization length £ also follows from the DMPK equation. 
It is given by £ = (f3N + 2 — /3)i, so it is larger than £ by a factor 
of order N. 
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FIG. 2 Sub-Poissonian shot noise in a disordered gold wire 
(dimensions 940 nm x 100 nm). At low currents the noise sat- 
urates at the level set by the temperature of 0.3 K. At higher 
currents the noise increases linearly with current, with a slope 
that is three times smaller tha n the va lue expected for Poisson 
statistics. Adapted from Ref. (|Hen99h . 



quantifies the deviation of the current fluctuations from a 
Poisson process (which would have F — 1). Since T <C 1, 
if all T„'s would be near T the current fluctuations would 
have Poisson statistics with F = 1. The bimodal dis- 
tribution (|2.14| instead gives sub-Poissonian shot noise 
(lBee92ft . 



F 



1 - 



JdTp(T)T 2 
J dTp(T)T 



(2.16) 



(The replacement of the sum over n by an integra- 
tion over T with weight p(T) is justified in the large- 
iV limit.) This one-third suppression of shot noise be- 
low the Poisso n value has been confirmed experimentally 
(|Ste96t lHen99h . see Fig. H 
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III. QUANTUM DOTS 

A. Level and wave function statistics 

Early applications of random matrix theory to con- 
densed matter physics were due to Gorkov and Eliashberg 
( Gor65 ) and to Denton, Miihlschlegel, and Scalapino 
(|Den7ir ). They took the Gaussian orthogonal ensem- 
ble to model the energy level statistics of small metal 
grains and used it to calculate quantum si ze effe cts on 
their thermodynamic properties. (See Ref. ( Hal86l ) for a 
review.) Theoretical justification came w ith the super- 
symmetric field theory of Efetov (Efc83), who derived 
the level correlation functions in an ensemble of disor- 
dered metal grains and showed that they agree with the 
RMT prediction up to an energy scale of the order of the 
inverse ergodic time h/r crg . 

Experimental evidence for RMT remained rare 
throughout the 1980's — basically because the energy 
resolution needed to probe spectral statistics on the scale 
of the level spacing was difficult to reach in metal grains. 
Two parallel advances in nanofabrication changed the sit- 
uation in the 1990's. 

One the one hand, it became possible to make electri- 
cal contact to in dividua l metal particles of diameters as 
small as 10 nm (|Del0lf ). Resonant tunneling through a 
single particle could probe the energy level spectr um with 
sufficient accuracy to test the RMT predictions ( Kue08[ ) 
(see Fig.©. 

On the other hand, semiconductor quantum dots be- 
came available. A quantum dot is a cavity of sub-micron 
dimensions, etched in a semiconducting two-dimensional 
electron gas. The electron wave length Xf — 50 nm at the 
Fermi energy in a quantum dot is two order of magnitudes 
greater than in a metal, and the correspondingly larger 
level spacing makes these systems ideal for the study of 
spectral statistics. The quantum dot may be disordered 
(mean free path I less than its linear dimension L) or it 
may be ballistic (I greater than L). RMT applies on en- 
ergy scales fi./T org ~ (hvp/L) min(l, l/L) irrespective of 
the ratio of I and L, provided that the classical dynamics 
is chaotic. 

Resonant tunneling through quantum dots has pro- 
vided detailed information on both the level spacing dis- 
tribution (through the spacing of the resonances) and on 
the wave functio n stati stics (through the peak height of 
the resonances ) (lAlhOOT i. For resonant tunneling through 
single-channel point contacts (tunnel probability T) the 
conductance peak height G max is related to the wave 
functio n intensities I±, 1% at the two point contacts by 
(|Bee9lh 



The Porter-Thomas distribution P(I) <x l^l 2 l e 



■PI/2 



G n 



TS hh 



h Ak B Th 



(3.1) 



(The intensities are normalized to unit average and 5 is 
the mean energy level spacing. The thermal energy kgT 
is assumed to be large compared to the width TS of the 
resonances but small compared to 6.) 



of (independently fluctuating) intensities I\ , I2 in the 
GOE {0 = 1) and GUE (0 = 2) then gives the peak 
height distribution (|Jal92HPri9"3h . 



P(a) = { ^)- 1/2 e-°., 
r(9) 1 g[K (g) + K^e-s, 



(3.2) 



with g = (8kBT/F5)(h/e 2 )G nia x and Bessel functions K , 
K\ . A comparison of this RMT prediction with available 
experimental data has shown a consistent agreement, 
with some deviations remaining that can be explained 
by finite-t empera ture effects and effects of exchange in- 
teraction (|Alh02t ). 



B. Scattering matrix ensembles 

In quantum dots, the most comprehensive test of RMT 
has been obtained by studying the statistics of the scat- 
tering matrix S rather than of the Hamiltonian H . The 
Hamiltonian H an d scattering m atrix S of a quantum 
dot are related by (|Bla9ll : lGot08[ ) 



S(E) = 1 - 2niW\E -H + inWW^^W 

_ t + inW^H - E)~ l W 
~ 1 - ittW(H- 



(3.3) 



The M x (Ni + N2) coupling matrix W (assumed to be 
independent of the energy E) couples the M energy lev- 
els in the quantum dot to N% + N2 scattering channels in 
a pair of point contacts that connect the quantum dot to 
electron reservoirs. The eigenvalue w n of the coupling- 
matrix product W^W is related to the transmission prob- 
ability r„ 6 [0, 1] of mode n through the point contact 
by 



An 2 w n M5 
{M5 + ir 2 w n )' 



(3.4) 



Eq. (|3.3[) is called the Wcidcnmullcr formula in the theory 
of chaotic scattering, because of pione ering work by Hans 
Wcidcnmullcr and his group (TMah69ri . 

A distribution function P(H) for the Hamiltonian H 
implies a distribution functional P[S(E)] for the scat- 
tering matrix S(E). For electrical conduction at low 
voltages and low temperatures, the energy may be fixed 
at the Fermi energy Ep and knowledge of the distribu- 
tion function P(Sq) of So = S(Ep) is sufficient. For the 
Hamiltonian we take the Gaussian ensemble, 



P(H) cx exp(-/5(7r/2<5) 2 M" 1 Tri/ 2 



(3.5) 



and we take the limit M — > 00 (at fixed 5, Ep, Y n ), 
appropriate for a quantum dot of size L Xp. The 
number of channels iVi , N2 in the two point contacts may 
be as small as 1, since the opening of the point contacts 
is typically of the same order as Xf. 
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FIG. 3 Data points: Integrated level spacing distribution of a single 10 nm diameter gold particle (barely visible in the 
micrograph as a white dot touching source and drain electrodes), measured by resonant tunneling in zero magnetic field and 
in a high magnetic field. The level spacings s are normalized by the mean level spacing 8 (equal to 0.23 meV in zero field 
and reduced to 0.12 meV in high fields due to splitting of the spin degenerate levels by the Zeeman effect). The measured 
distributions are compared with Wigner's RMT prediction: P(s) oc s 13 exp(— cgs 2 ) (with ci = 7r/4, C2 = 4/-7T, C4 = 64/97r), for 
the Gaussian orthogonal ensemble (GOE, /3 = 1), the Gaussian unitary ensemble (GUE f3 = 2) and the Gaussian symplectic 
ensemble (GSE, = 4). The Poisson distribution P(s) oc e~ s of uncorrelated levels is also shown. A magnetic field causes 
a transition from the symplectic ensemble in zero field (preserved time reversal symmetry, broken spin rotation symmetry 
due to the strong spin-orbit coup ling in gold), to the unitary ensemble in high fields (broken time reversal and spin rotation 
symmetries). Adapted from Rcf. (Kuc08). 



As derived by Brouwer ( Bro95f) . Eqs. f|3.3|) and (|3.5p 
together imply, in the large- M limit, for So a distribution 
of the form 

P(S ) oc |Det(l - S^So)]-* 31 * 1 - 01 * 2 - 2 ^, (3.6) 

known as the Poisson kernel (jHuafi3l : iLewQll : lDor92f L 

The average scattering matrix 4 S = j dSo SqP(Sq) in 
the Poisson kernel is given by 



MS 



MS + tt 2 WW 



(3.7) 



The case of ideal coupling (all r„ 's equal to unity) is of 
particular interest, since it applies to the experimentally 
relevant case of ballistic point contacts (no tunnel barrier 
separating the quantum dot from the electron reservoirs) . 
In view of Eq. (|3.4[) one then has S = 0, hence 



P(S ) 



constant. 



(3.8) 



This is the distribution of Dyson's circular ensemble 
(1DVS62TI . irst a pplied to quantum scattering by Blumcl 
and Smilansky ( Blu90l ). 

The circular ensemble of scattering matrices implies 
for the min( A^i , N2) nonzero transmission eigenvalues the 
distribution (|Bee97l ) 



pirn)* n iT n -r m |^n r , 



1/3(1^2-^11 + 1-2/^) 



n<m 



(3.9) 



4 The average S is defined by integration over the unitary group 
with Haar measure dSo, unconstrained for /3 = 2 and subject 
to the constraints of time reversal symmetry for (3 = 1 (when 
S is symmetric) or symplectic symmetry for f3 = 4 (when S is 
self-dual). For m ore information on integration over the unitary 
group, see Refs. (iBcc97; Guh9Sh. 



This distribution is of the form (|2.3p . with the logarith- 
mic repulsion (|2.4[) . There are no nonlogarithmic correc- 
tions in a quantum dot, unlike in a quantum wire. 



C. Conductance distribution 

The complete probability distribution of the conduc- 
tance G = G Q J2n=} Nl ' N2)T n follows directly from Eq. 
(|3.9[) in the cas e N\ = N2 = 1 of single-channel ballistic 
point contacts 



P(G) = 2^(G/G ) 



-1+/3/2 



< G < G . 



(3.10) 



This strongly non-Gaussian distribution rapidly ap- 
proaches a Gaussian with increasing Ni = N2 = N. Ex- 
periments typically find a conductance distribution which 
is close r to a Gaussian even in the single-channel case 
(Hui98), due to thermal averaging and loss of phase co- 
herence at finite temperatures. 

In the limit N — » 00 the variance of the Gaussian is 
given by the RMT result IpTf]) for UCF — without any 
corrections since the eigenvalue repulsion in a quantum 
dot is strictly logarithmic. The experiment value in Fig. 
|4] is smaller than this zero-temperature result, but the 
factor-of-two reduction upon application of a magnetic 
field ((3 = 1 — > P = 2) is quite well preserved. 

Without phase coherence the conductance would have 
average GqN/2, corresponding to two A^-mode point con- 
tacts in series. Quantum interference corrects that aver- 
age, (G) = G0A72 + 5G. The correction SG in the limit 
A^ — > 00, following from the circular ensemble, equals 



M-i(l-|'| C „. 



(3.11) 
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FIG. 4 Variance of the conductance of a quantum dot at 
30mK, as a function of magnetic field. The inset shows 
an electron micrograph of the device, fabricated in the two- 
dimensional electron gas of a GaAs/AlGaAs heterostructure. 
The black rectangle at the center of the inset is the quantum 
dot, the gray regions are the gate electrodes on top of the 
heterostructure. Electrons can enter and exit the quantum 
dot through point contacts at the top and right corner of the 
rectangle. The side of the rectangle between these two corners 
is distorted to generate conducta nce fluc tuations and obtain 
the variance. Adapted from Rcf. (Cha95). 
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FIG. 6 Sub-Poissonian shot noise in a quantum dot at 
270 mK. The slope at high currents corresponds to a one- 
quarter Fano factor, as predicted by RMT. The inset shows 
an electron micrograph of the device. The quantum dot is 
contained between point contacts A and B. (The gate labeled 
C is no t operative in this experiment.) Adapted from Ref. 
(|Obe0lt) . 




FIG. 5 Magnetoconductance at 50 mK, averaged over 48 
quantum dots. The minimum around zero magnetic field 
is the weak localization effect. The inset shows the geom- 
etry of the quantum dots, which are fabricated in the two- 
dimensional electro n gas of a GaAs/AlGaAs heterostructure. 
Adapted from Ref. (|Cha94l ). 



The quantum correction vanishes in the presence of a 
time-reversal-symmctry breaking magnetic field (/? = 2), 
while in zero magnetic field the correction can be nega- 
tive (j3 = 1) or positive (j3 — 4) depending on whether 
spin-rotation-symmetry is preserved or not. The nega- 
tive quantum correction is called weak localization and 
the positive quantum correction is called weak antilocal- 



ization. An experimental demonstration ( Cha94l ) of the 
suppression of weak localization by a magnetic field is 
shown in Fig. [5] The measured magnitude 6G of the peak 
around zero magnetic field is 0.2 Go, somewhat smaller 
than the fully phase-coherent value of 4 Go ■ 



D. Sub-Poissonian shot noise 

For Ni — A r 2 = N ^> 1 the density of transmission 
eigenvalues for a quantum dot, following from Eq. (|3.9[) . 
has the form 



P(T) = 



N 



1 



7T Vt^/T^t 



(3.12) 



It is different from the result (|2. 14(1 for a wire, but it has 
the same bimodal structure: While the average transmis- 
sion T = 1/2, the eigenvalue density is peaked at zero and 
unit transmission. 

This bimodal structure can be detected as sub- 
Poissonian shot noise. Instead of Eq. (|2.16p one now has 
(|Jal94D 



F 



JdTp(T)T 2 



1 



/ dTp(T)T 4' 
An experimental demonstration is shown in Fig. [5] 



(3.13) 
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FIG. 7 Thermopower distribution for /3 = 1 (\B\ < 40mT) and (3 = 2 (\B\ > 50 mT). Experimental results (dots), RMT 
results (solid line), and Gaussian fit (dashed line) are compared. The inset shows the experimental layout. The crosses denote 
Ohmic contacts to the two-dimensional electron gas and the shaded areas denote gate electrodes. The heating current is applied 
between I\ and I2, while th e thermovoltage is measured between Vi and Vi. The quantum dot is indicated by the dotted area. 
Adapted from Ref. (|God99t ). 



E. Thermopower distribution 

Knowledge of the distribution of the scattering matrix 
S(E) at a single energy E = Ep is sufficient to determine 
the conductance distribution, but other transport prop- 
erties require also information on the energy dependence 
of S. The thermopower V (giving the voltage produced 
by a temperature difference at zero electrical current) is a 
notable example. Since V oc dlnG/dE, we need to know 
the joint distribution of S and dS/dE at Ep to determine 
the distribution of V . 

This problem goes back to the early days of RMT 
(IWig55l : ISmi60l ). in connection with the question: What 



is the time delay experienced by a scattered wave packet? 
The delay times r n are the eigenvalues of the Hermi- 
tian matrix product Qws = —ifiS^dS/dE, known as 
the Wigner-Smith matrix in the context of RMT. [For 
applica tions in other contexts, see Refs. ljDas69l ; lBla9lt 
lGot08f ).] The solution to the problem of the joint dis- 
tribution of S and dS /dE (fo r S in the circular ensem- 
ble) was given in Ref. (|Bro97l ). The symmetrized matrix 
product 



Q = -ihs-v^s- 1 ' 2 

dE 



(3.14) 



has the same eigenvalues as Qws> but unlike Qws wa s 
found to be statistically independent of S. The eigenval- 
ues of Q have distribution 



P(bn})* I]>-^f Il7), 

i<j k 

-in = l/r n > 0. 



P TT ^( Ar i + Ar 2)/2 e - / 3r„ 7fc /2 



(3.15) 



The Heisenberg time tjj = 2-kTi/5 is inversely propor- 
tional to the mean level spacing 5 in the quantum dot. 
Eq. (|3 . 1 5|) is known in RMT as the Laguerre ensemble. 

The thermopower distribution following from the La- 
guerre ensemble is strongly non-Gaussian for small N± = 



A 2 = N. For N = 1 it has a cusp at V — when 
(3=1 and algebraically decaying tails oc |7 : '|~ 1 ~ /3 In \V\. 
Significant deviatio ns from a Gaussian are seen in the 
experiment (|God99H shown in Fig. [71 for N = 2. 



F. Quantum-to-classical transition 

RMT is a quantum mechanical theory which breaks 
down in the classical limit h — > 0. For electrical con- 
duction through a quantum dot, the parameter which 
governs the quantum-to-classical transition is the ratio 
TE/Vdwcii of Ehrcnfest time and dwell time 

The dwell time Td wo ii is the average time an electron 
spends inside the quantum dot between entrance and exit 
through one of the two A^-mode point contacts. It is given 
by 



Tdwoii = irh/N5. 



(3.16) 



The Ehrenfest time Tp is the time up to which a wave 
packet follows classical e quations of m otion, in accord 
with Ehrenfest 's theorem ( Ber78l : IChi7lh . For chao tic dy - 
namics with Lyapunov exponent 5 a, it is given by (|Sil03h 

t e = a- 1 max[0, \n(NW/A 1/2 )} . (3.17) 

Here A is the area of the quantum dot and W the width 
of the A'-mode point contacts. 

The RMT result F = 1/4 holds if t e < r dweU . For 
longer T p, the Fano factor is suppressed exponentially 
(|Aga00h . 



F 



4 e 



(3.18) 



The Lyapunov exponent a of chaotic motion quantifies the ex- 
ponential divergence of two trajectories displaced by a distance 
Ax(t) at time t, according to Ax(t) = Ai(0)e al . 
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FIG. 8 Dependence of the Fano factor F of a ballistic chaotic 
quantum dot on the average time that an electron dwells in- 
side. The data points with error bars are measured in a quan- 
tum dot, the solid curve is the theoretical prediction (13.18[) 
for the quantum-to-classical transition (with Ehrenfest time 
te = 0.27ns as a fit parameter). The dwell time (|3.16|l is 
varied experimentally by changing the number of modes N 
transmitted through each of the point contacts. The inset 
shows graphically the sensitivity to ini tial con ditions of the 
chaotic dynamics. Ada pted fro m Ref. (|Aga00l) . with experi- 



mental data from Ref. l|Obe02l ) 



This equation expresses the fact that the fraction 1 — 
e -TE/T dwoll Q £ e i ec t rons that stay inside the quantum dot 
for times shorter than te follow a deterministic classical 
motion that does not contribute to the shot noise. RMT 
applies effectively only to the fraction e _TE / Tdwdl of elec- 
trons that stay inside for times longer than te- The shot 
noise suppression (|3. T8[) is plott ed in Fi g. El together with 
supporting experimental data ( Obe02l ). 



equation 



Det [1 - a(E) 2 S(E)S(-E)*] = 0, 



a(E) = | 



E 2 



(4.1) 



The scattering matrix S (at an energy E measured rel- 
ative to the Fermi level) describes the coupling of the 
quantum dot to the superconductor via an iV-mode point 
contact and is related to the Hamiltonian H of the iso- 
lated quantum dot by Eq. (|3.3| . At low energies E <C A 
the energy levels can be obtained as the eigenvalues Ei 
of the effective Hamiltonian 



H 

-nW*W^ 



-nWW T 
-H* 



(4.2) 



The Hcrmitian matrix H e f{ is antisymmetric under the 
combined ope ration of charge conjugation (C) and time 
inversion (T) {Alt96|): 



-a y H* B a y , 



-i 

1 



(4.3) 



(An M x M unit matrix in each of the four blocks of a y is 
implicit.) The CT-antisymmetry ensures that the eigen- 
values lie symmetrically around E = 0. Only the positive 
eigenvalues are retained in the excitation spectrum, but 
the presence of the negative eigenvalues is felt as a level 
repulsion near E = 0. 

As illustrated in Fig. [5] (left panel), the unique feature 
of the proximity effect is that this level repulsion can 
extend over energy scales much larger than the mean 
level spacing S in the isolated quantum dot — at least 
if time reversal symmetry is not broken. A calculation 
of the density of states (p(E)) = (£\ S(E - Ei)) of H eS , 
averaged over H in the GOE, produces a square root 
singularity in the large- N limit: 



(p(E)) - - 

7T 



1 E-E, 



gap 



A 3 

gap 



E 



E N - 



oo, (4.4) 



IV. SUPERCONDUCTORS 



A. Proximity effect 



Fig. [9] (lower right panel) shows a quantum dot with 
superconducting electrodes. Without the superconduc- 
tor the energy spectrum of an ensemble of such quantum 
dots has GOE statistics. The proximity of a supercon- 
ductor has a drastic effect on the energy spectrum, by 
opening up a gap at the Fermi level. The RM T of this 
proxim ity effect was developed in Ref. (|Mel96t ) ( see Ref. 
(|Bee05h for a review). 

A quantum dot coupled to a superconductor has a dis- 
crete spectrum for energies below the gap A of the su- 
perconductor, given by the roots of the determinantal 



If the point contact between quantum dot and supercon- 
ductor is ballistic (T n = 1 for n = 1, 2, . . . N) the two 
energies E gap and A gap are given by (|Mer96T) 



2n 



= 0.048^5, A gap = 0.068 N 1/3 S. 

(4.5) 

(Here 7 = \(\/Z — 1) is the golden number.) The gap 
-Egap in the spectrum of the quantum dot is larger than 
8 by factor of order N. 



B. Gap fluctuations 



The value (|4. 5|) of the excitation gap is representative 
for an ensemble of quantum dots, but each member of the 
ensemble will have a smallest excitation energy E\ that 
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FIG. 9 Left panel: Average density of states (scaled by the Thouless energy Et = NS/Att) of a quantum dot coupled by a 
ballistic iV-mode point contact to a superconductor. The histogram is a numerical calculation for the ge ometry indicated in 
the inset (with N = 20), while the curve is the analytical prediction from RMT. Adapted from Ref. (JMej96j). Upper rig ht pane l: 
Probability distribution of the lowest excitation energy Ex, rescaled as x = {Ex - £ g ap)/A gap . Adapted from Ref. (|VavOj) . 
Lower right panel: Quantum dot (central square of dimensions 500 nm x 500 nm) fabricated in an InAs/AlSb heterostructure 
and con tacted by four superconducting Nb electrodes. Device made by A.T. Filip, Groningen University. Figure from Ref. 
(|Bee05l V 



will be slightly different from E gap . The distribution of is governed by the distribution 



the gap fluctuations is ide ntical upon rescaling to the 
known dis tribution (lTra94) of t he lowest eigenvalue in 
the GOE (|Vav0ll ; [OstOll iLamOll ). Rescaling amounts to 



a change of variables from E\ to x = {E\ — _E gap )/A g 
where -E gap and A gap parameterize the square-root de- 



pendence (|4.4p . The probability distribution P{x) of the 
rescaled gap fluctuations is shown in Fig. [9] (upper right 
panel). The gap fluctuations are a mesoscopic, rather 
than a microscopic effect, because the typical magnitude 
A gap ~ -E g ap<5 2 / 3 of the fluctuations is 8 for E gap ^> 8. 
Still, the fluctuations are small on the scale of the gap 
itself. 



P({E n })^H\Ef-E^l[\E k fe- 

i<j k 



(4.6) 



related to the Lagucrrc ensemble by a change of variables 
{E 2 — » x n ). (The coefficient c is fixed by the mean level 
spacing in the isolated quantum dot.) The density of 
states near zero energy vanishes as \E\°. Two more cases 
are possible when spin-rotation symmetry is broken, so 
that in total the three Wigner-Dyson symmetry classes 
without superconductivity are expanded to four symme- 
try classes as a consequence of the CT-antisymmetry. 



D. Quantum-to-classical transition 



C. From mesoscopic to microscopic gap 

The mesoscopic excitation gap of order N8 induced by 
the proximity to a superconductor is strongly reduced 
if time reversal symmetry is broken by application of a 
magnetic field ([3 = 2). Because the repulsion of levels 
at ±E persists, as demanded by the CT-antisymmetry 
(|4.3p . a microscopic gap around zero energy of order 8 
remains. An alternative way to reduce the gap from N8 
to 8, without breaking time reversal symmetry {(3 = 1), 
is by contacting the quantum dot to a pair of supercon- 
ductors with a phase difference of n in the o rder pa ram- 
eter. As shown by Altland and Zirnbauer ( Alt96l ) . the 
level statistics near the Fermi energy in these two cases 



The RMT of the proximity effect discussed so far 
breaks down when the dwell time (13.16)) becomes shorter 
than the Ehrenfest time (|3.17p (Lod98) . In order of mag- 
nitude, 6 the gap equals -E gap ~ min{Ti/TE,h/Tdwen)- In 
the classic al limit te — > oo, the density of states is given 
by (|Sch99h 



= 2 {ttE t /E) 2 cosh{irE T /E) 
{P( 8 sinh 2 {ttEt/E) 



(4.7) 



More precisely, the gap crosses over between the RMT limit (14.51 



for te < S Td we n and t he limit E g 
(Vav03; Bcc05; Kui09). 



ttTi/2te for t e > T dwo n 
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FIG. 10 Comparison of the density of states (|4.7[1 with the 
RMT result (|4.4[) . These are the two limiting results when 
the Ehrenfest time te is, respectively, much larger or mu ch 
smaller than the mean dwell time Tdweii- From Ref. (|Bee05l ). 



This is an example of an application of RMT to clas- 
sical optics, because the spectral statistics of a cavity is 
determined by the Maxwell equations of a classical elec- 
tromagnetic wave. (More applications of this type, i n- 
cluding also sound waves, are reviewed in Ref. ()Kuh05lU 
An altogether different type of application of RMT ap- 
pears in quantum optics, when the photon and its quan- 
tum statistics play an essential role. Selected applica- 
tions of RMT to both classical and quantum optics are 
presented in the following sections. The emphasis is on 
topics that do not have an immediate analogue in elec- 
tronics, either because they cannot readily be measured 
in the solid state or because they involve aspects (such 
as absorption, amplification or bosonic statistics) that do 
not apply to electrons. 

Some of the concepts used in this Chapter were intro- 
duced in the previous Chapter on applications of RMT 
to condensed matter physics, in particular in Sees. III. Bt 
HTUl andQTTBj 



with Et = NS/Att the Thouless energy. The density of 
states (|4.7I) (plotted in Fig. [10]) is suppressed exponen- 
tially oc e ~* E T/E at the F erm i i eve i (£ _> q), but there 

is no gap. 

To understand the absence of a true excitation gap in 
the limit te — > oo, we note that in this limit a wave 
packet follows a classical trajectory in the quantum dot. 
The duration t of this trajectory, from one reflection at 
the superconductor to the next, is related to the energy 
E of the wave packet by E ~ h/t. Since t can become ar- 
bitrarily large (albeit with an exponentially small prob- 
ability e _t/Td ™ 011 ), the energy E can become arbitrarily 
small and there is no gap. 



Part II 

Classical and Quantum 
Optics 

V. INTRODUCTION 

Optical applications of random matrix theory came 
later than electronic applications, perhaps because ran- 
domness is much more easily avoided in optics than it is 
in electronics. The variety of optical systems to which 
RMT can be app l ied inc reased substantially with the re- 
alization ||Boh84 lBer85t ) that randomness is not needed 
at all for GOE statistics of the spectrum. Chaotic dy- 
namics is sufficient, and this is a generic property of res- 
onators formed by a combination of convex and concave 
surface elements. As an example, we show in Fig. [11] 
the Wigner level spacing distribution measured in a mi- 
crowave cavity with a chaotic shape. 



VI. CLASSICAL OPTICS 

A. Optical speckle and coherent backscattering 

Optical speckle, shown in Fig. 1121 is the random inter- 
ference pattern that is observed when coherent radiation 
is transmitted or reflected by a random medium. It has 
been much studied since the discovery of the laser, be- 
cause the speckle pattern carries information both on the 
coherence properties of the radi ation an d on microscopic 
details of the scattering object ( Goo07l ). The superposi- 
tion of partial waves with randomly varying phase and 
amplitude produces a wide distribution P(I) of intensi- 
ties I around the average /. For full coherence and com- 
plete randomization the distribution has the exponential 
form 



P(I) = J" 1 exp(-J/J), I > 0. 



(6.1) 



For a descr iption of speckle in the framework of 
RMT (Mel88a), it is convenient to enclose the scatter- 
ing medium in a wave guide containing a large number 
N of propagating modes. The reflection matrix r is then 
aniVxJV matrix with random elements. Time-reversal 
symmetry (reciprocity) dictates that r is symmetric. De- 
viations of r from unitarity can be ignored if the mean 
free path I is much smaller than both the length L of the 
scattering medium and the absorption length l a . The 
RMT assumption is that r is distributed according to 
the circular orthogonal ensemble (COE), which means 
that r = UU T with U uniformly distributed in the group 
U(N) of N x N unitary matrices. 

In this description, the reflected intensity in mode n 
for a wave incident in mode m is given by I nm — \r nm \ . 
The intensity distribution can be easily calculated in the 
limit N — > oo, when the complex matrix elements r nm 
with n < m have independent Gaussian distributions of 
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FIG. 11 Histogram: distribution of spacings s of eigenfrequencies measured in the chaotic microwave resonator shown at the 
right. (The resonator has superconducting walls, to minimize absorption.) The spacing distribution is close to the Wigner 
distribution P(s) oc s exp(— 7rs 2 /45 2 ) [solid line] of the GOE, and far from the Poisson distribution P(s) oc e _S//<5 [dashed line] 
of uncorrelated eigenfrequencies. The mean spacing has been set to 8 = 1, a nd non-chaotic "bouncing-ball" resonances have 
been eliminated from the experimental histogram. Adapted from Ref. (|Alt97h . 





intensity (arb. units) 



FIG. 12 Left panel: Speckle pattern produced by a laser beam behind a diffusor (full scale 45mrad x 45mrad). The vertical 
bar indicates the color coding of the intensity, in arbitrary units. The average angular opening angle Sa ~ 1.3 mrad of a bright 
or dark spot (a "speckle") is equal to \/ttR, with A = 830 nm the wave length and 7? = 200 fj,m the radius of the illuminated 
area on the diffusor. The envelope of the intensity pattern reflects the 18 mrad opening angle of the directional scattering from 
this type of diffusor. The intensity distribution P(I) of the speckle pattern measured inside the white square is plotted in the 
right panel, and compared with the exponential distribution (|6.1[) (straight line in the semi- logarithmic plot). Figure courtesy 
of M.P. van Exter. 



zero mean and variance 7 



(l»w| a ) 



JV 

du ]T u nk u mk ,u I 

"W fe,fe'=i 



1 



N + 1 



(6.2) 

The resulting distribution of I nm in the large-iV limit 
has the exponential form ()6. 1|) . with an average inten- 
sity I nm — (1 + S nm )N which is twice as large when 
n = m than when n ^ m. This doubling of the average 
reflected intensity at t he angl e of incidence is the coherent 
backscattering effect (]Akk07| ) , illustrated in Fig. [T51 



7 For an intro duction to such integrals over the unitary group, see 
Ref. jBee97T) . The factor N + 1 in the denominator ensures that 
5Zm=i l r n>n| 2 = 1) as required by unitarity, but the difference 
between N and N + 1 can be neglected in the large- ./V limit. 



The RMT assumption of a COE distribution of the 
reflection matrix correctly reproduces the height of the 
coherent backsca t tering p eak, but it cannot reproduce 
its width (|Akk88l : lMel88bt ). The Kronecker delta in Eq. 
(16. 2|) would imply an angular opening Sa ~ 1/kW of the 
peak (for light of wave number A: in a wave guide of width 
W). This is only correct if the mean free path I is larger 
than W. In a typical experiment I -C W and the angular 
opening is Sa ~ 1/kl (as it is in Fig.[Tr 



B. Reflection from an absorbing random medium 

An absorbing medium has a dielectric constant e with a 
positive imaginary part. The intensity of radiation which 
has propagated without scattering over a distance L is 
then multiplied by a factor e~ aL . The decay rate a > 
at wave number k is related to the dielectric constant by 
a = 2k Im y/e. 
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FIG. 13 Measurement of coherent backscattering from a ZnO 
powder. The sample is rotated to average the reflected inten- 
sity, which is plotted against the scattering angle. The mea- 
sured peak due to coherent backscattering is superimposed on 
the diffuse scattering intensity (dashed curve, normalized to 
unity in the backscattering direction at zero angle). The rel- 
ative height of the peak is a factor-of-two, in accord with Eq. 
(|6.2|) . The angular width is of order 1/kl ~ 40 mrad, for wave 
length A = 2ir/k = 514 nm and mean free path I = 1.89 /im. 
The inset shows the optical setup. Use of a beam splitter per- 
mits detection in the backscattering direction, which would 
otherwise be bloc ked by the incident laser beam. Adapted 
from Ref. (|Wie95h . 



The absence of a conservation law in an absorb- 
ing medium breaks the unitarity of the scattering ma- 
trix. The circular orthogonal ensemble, of uniformly dis- 
tributed symmetric unitary matrices, should therefore be 
replaced by another ense mble. The app ropriate ensem- 
ble was derived in Refs. ( Bee9 6; Bru96), for the case of 
reflection from an infinitely long absorbing wave guide. 
The result is that the N eigenvalues R n E [0,1] of the 
reflection matrix product rr^ are distributed according 
to the Laguerre orthogonal ensemble, after a change of 
variables to A„ = R n (l — i? rl ) _1 > 0: 



P({K}) oc Yl \Xj - Ai| l[exp[-al(N + l)X k ], (6.3) 



The distribution (|6.3p is obtained by including an ab- 
sorption term into the DMPK equation (l2T2p. This loss- 
drift-diffusion equation has the form (|Bee96l ; lBru96l ) 



JY 



d 



,dP 



dL (3N + 2 - ^ d\ n 

n—1 



jJ!-^+al(/3N + 2-l3)P 



d\ n J 

with J = X\\\j - M 13 



(6.4) 



i<3 



The drift-diffusion equation (|2.12p considered in the elec- 
tronic context is obtained by setting a — and trans- 
forming to the variables x n = sinh 2 X n . 



In the limit L — > oo we may equate the left-hand-side 
of Eq. (|6.4[) to zero, and we arrive at the solution (|6.3p 
for [3 — 1 (unbroken time reversal symmetry). More 
generally, for any (3, the distribution of the R n 's in the 
limit L — > oo can be written in the form of a Gibbs 
distribution at a fictitious temperature /3 , 



P({R n }) oc exp [-P(J2 <R t , Rj)+Y, V< y R i) 



u{R,R') = - \n\R-R'\, 
2 



V(R) 



N-1 + 



alR 
l-R 



ln(l - R) 



(6.5) 



(6.6) 



The eigenvalue interaction potential u(R, R') is logarith- 
mic. This can be contrasted with the nonlogarithmic 
interaction potential in the absence of absorption, dis- 
cussed in Sec. III.B1 Because R n = 1 — T n without ab- 
sorption, the interaction potential (|2.9|) of that section 
can be written as 



u(R, R') = 



■\\n\R- 



■R'\-l\n\ 



: |, with R = tanh x. 

(6.7) 

As calculated in Ref. ( Mis96f ). the change in interaction 
potential has an observable consequence in the sample- 
to-sample fluctuations of the reflectance 



N 



1Z = Tr rr 



22 Rr, 



(6.8) 



With increasing length L of the absorbing disordered 
waveguide, the variance of the reflectance drops from the 
value VarT?. = 2/15/3 associated with the nonlogarithmic 
interaction (|6~71) [cf. Eq. fljlj ], to the value Varft = 1/8/3 
for a logarithmic interaction [cf. Eq. (|2.7p ]. The crossover 
occurs when L becomes longer than the absorption length 
la = \/l/o~a, in the large- regime A^ » \j\fal 1. 



C. Long-range wave function correlations 

The statistics of wave function intensities / = | \I/ (t*) | 2 
in a chaotic cavity is described by the Porter-Thomas 
distribution (|Por65tl . 



P(I) = (2^J)- 1 / 2 /- 1 / 2 cxp(-//27), J>0, (6.9) 

with / the average intensity. Eq. (|6.9p assumes time re- 
versal symmetry, so ^ is real (symmetry index (3 = 1). 
An experimental demonstration in a microwave resonator 
is shown in Fig. [I4j 

In the context of RMT, the distribution (|6.9p follows 
from the GOE ensemble of the real symmetric M x M 
matrix H (the effective Hamiltonian) , which determines 
the eigenstates of the cavity. The intensity / corresponds 
to the square of a matrix element O nm of the orthogonal 
matrix which diagonalizes H, where the index n labels 
a point in discretized space and the index m labels a 
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FIG. 14 Comparison of the Porter-Thomas distribution (16. 9[) 
[dashed line] of wave function intensities / in the GOE, with 
the intensity distribution measured on the two-dimensional 
microwave cavity shown in the inset. (The average intensity 
has been set to / = 1.) Adapted from Ref. (|Kud95l ). 



particular eigenstate. In the large-M limit the matrix 
elements of O have a Gaussian distribution, which implies 
Eq. for the distribution of / = 2 m . 

Different matrix elements O nm and O n i m are indepen- 
dent, so the wave function has no spatial correlations 
in the RMT description. This is an approximation, but 
since the act ual co rrelations decay on the scale of the 
wave length ( Ber77l ). it is accurate to say that there are 
no long-range wave function correlations in a chaotic cav- 
ity. 

The same absence of long-range correlations applies if 
time reversal symmetry is fully broken, by the introduc- 
tion of a s ufficient ly strong magneto-optical element in 
the cavity ( Sto99l ). The intensity distribution changes 
from the Porter-Thomas distribution (|6.9p to the expo- 
nential distribution (|6.ip . but spatial correlations still 
decay on the scale of the wave length. Partially bro- 
ken time reversal symmetry, however, has the striking ef- 
fect of introducing wave function correlations that persist 
throughout the entire cavity. T his wa s discovered theo- 
retically by Fal'ko and Efetov (|Fal94l ) for the crossover 
from GOE to GUE. 

An altogether different way to partially break time re- 
versal symmetry is to open up the cavity by attaching a 
pair of A^-mode leads to it, an d to ex cite a traveli ng wave 
from one lead to the other ( Pni96f ). Brouwer ( Bro03f ) 
found that, if A" is of order unity, the traveling wave 
produces relatively large long-range wave function corre- 
lations inside the cavity. As shown in Fig. [T5J these cor- 
relations have been measured in a microwave resonator 
(|Kim05t ). 

Partially broken time reversal symmetry means that 
a wave functio n ^(r) is neither real nor fully complex. 
Following Ref. (lLan97T) . the crossover from real to fully 
complex wave functions is quantified by the phase rigidity 



(6.10) 



p Jdr\^{r)\ 2 ' 
A real wave function has p = 1 while a fully complex 



wave function has p — 0. 

As \p\ decreases from 1 to 0, the intensity distribution 
crosses over from the Porter-Thomas distribution (16.91) to 
the exponential distribution (|6.1[) . according to (|Pni96T ) 



P(I\P) 



1 



■. exp 



I/I 
l-\p\ 



W/i 



1 



\p\ 2 , 

(6.11) 

(The function Iq is a Bessel function.) The notation 
P{I\p) indicates that this is the intensity distribution for 
an eigenstate with a given value of p. The distribution 
P(p) of p among different eigenstates, calculated in Ref. 
( Bro03l ) , is broad for N of order unity. 

For any given phase rigidity the joint distribution of 
the intensities / = I(r) and I' = I(r') factorizes if 
k\r — r'\ 1. The long-range correlations appear upon 
averaging over the broad distribution of phase rigidities, 
since 



P(I,I') = J d P P(p)P(I\p)P(l'\p) 



(6.12) 



no longer factorizes. 



D. Open transmission channels 



The bimodal transmission distr ibution (12.14p , first ob- 
tained by Dorokhov in 1984 ( Dor84h . tells us that a frac- 
tion l/L of the transmission eigenvalues through a ran- 
dom medium is of order unity, the remainder being ex- 
ponentially small. A physical consequence of these open 
channels, discussed in Sec. III. CI is the sub-Poissonian 



shot noise of electrical current ( Bee92t ). As expressed by 
Eq. ()2.16|) . the shot noise power is reduced by a fac- 
tor 1 — 2/3 = 1/3, because the spectral average T 2 of 
the transmission eigenvalues is 2/3 of the average trans- 
mission T = l/L. If all transmission eigenvalues would 
have been close to their average, one would have found 
T 2 /T ~ l/L <C 1 and the shot noise would have been 
Poissonian. 

The observation of sub-Poissonian shot noise is evi- 
dence for the existence of open transmission channels, but 
it is indirect evidence — because a theory is required to 
interpret the observed shot noise in terms of the transmis- 
sion eigenvalues. In fact, one can alternatively interpret 
the sub-Poissonian shot noise in terms of a semiclassical 
theory t hat ma kes no reference at all to the transmission 
matrix ( Nag92l ). 

A direct measurement of the ratio T 2 /T would require 
the preparation of a specific scattering state, which is 
not feasible in electronics. In optics, however, this is a 
feasible experiment — as de monstrated very recently by 
Vellekoop and Mosk (|yel08). By adjusting the relative 
amplitude and phase of a superposition of plane waves, 
they produced an incident wave with amplitude -E™ = 
t* mon in mode n = 1,2,... N (for N ~ 10 4 ). The index 
mo corresponds to an arbitrarily chosen "target speckle" 
behind a diffusor, located at the center of the square 
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FIG. 15 Left panel: Distribution of the intensity |^(r)| 2 of a traveling wave at a fixed frequency in the open two-dimensional 
chaotic microwave cavity shown in the inset (dimensions 21 cm x 18 cm). The wave function $ is the component of the electric 
field perpendicular to the cavity (normalized to unit average intensity), for a wave traveling from the right to the left lead. The 
measured values (histogram) are compared with the distribution (16. lip (solid curve), fitted to a phase rigidity \p\ 2 = 0.5202. 
The grey scale plot in the inset shows the spatial intensity variations, with black corresponding to maximal intensity. Right 
panel: Correlator of squared intensity, for a single mode in both the right and left leads. The histogram shows the measured 
correlat or, averaged over position in the cavity and frequency of the traveling wave. Th e solid cu rve is the theoretical prediction 
(|Bro03T ), which tends to the nonzero limit 0.078 for k\r - r'\ > 1. Adapted from Ref. (|Kim05h . 




FIG. 16 Right panels: Speckle pattern (area 30 /J,m x 30 (im) behind a diffusor (a 11.3 /im layer of ZnO particles with mean 
free path I — 0.85 /im), for a random incident wave front (top) and for a wave front optimized to couple to open transmission 
channels (bottom). The intensity of the bright speckle at the center in the bottom panel is a factor of 750 greater than the 
background. Left panel: Intensity profile, integrated over the y-direction to average out the speck le pattern. The optimized 
wave front (red) has a peak, which the random wave front (black) lacks. Adapted from Ref. (|Vel08h . 



speckle pattern in Fig. 1161 The transmitted wave has 
amplitude 



intensity is a factor of order Nl/L ~ 10 3 , as observed. 
The total transmitted intensity is 



E° 



7- P 11 



(ttt) 



rnm.Q ■ 



(6.13) 



(6.14) 



As shown in Ref. ( Vel08l ). this optimized incident wave 
front can be constructed "by trial and error" without 
prior knowledge of the transmission matrix, because it 
maximizes the transmitted intensity at the target speckle 
(for a fixed incident intensity). The optimal increase in 



The average transmitted intensity, averaged over the tar- 
get speckle, gives the spectral average T 2 , 

7^ = — V 7 out = — Tr (t^ f = T 2 ". (6.15) 

mo 
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The average incident intensity is simply I m = 
N~ lr Trtt' = T, so the ratio of transmitted and incident 
intensities gives the required ratio of spectral averages, 
J out // m = T 2 /T. The experimental results are consis- 
tent with the value 2/3 for this ratio, in accord with the 
bimodal transmission distribution (|2.14|) . 



VII. QUANTUM OPTICS 

A. Grey-body radiation 

The emission of photons by matter in thermal equilib- 
rium is not a series of independent eve nts. Th e text- 
book example is black-body radiation (|Man95l ): Con- 
sider a system in thermal equilibrium (temperature T) 
that fully absorbs any incident radiation in N propagat- 
ing modes within a frequency interval Slo around w. A 
photodetector counts the emission of n photons in this 
frequency interval during a long time t ^> 1/Slu. The 
probability distribution P(n) is given by the negative- 
binomial distribution with v — Ntduj /2tt degrees of free- 
dom, 



The relation between F(£) and S is 



P(n) oc 



n + v — 1 



cxp ( — n huj J &b T ) • 



(7.1) 



The binomial coefficient counts the number of parti- 
tions of n bosons among v states. The mean photocount 
n = vf is proportional to the Bose-Einstein function 



/(w,T)=[eoq>(7iw/AtoT)-l]- 



(7.2) 



In the limit n/v — > 0, Eq. (|7.1|) approaches the Poisson 
distribution P(n) oc n n /n\ of independent photocounts. 
The Poisson distribution has variance Var n — n equal to 
its mean. The negative-binomial distribution describes 
photocounts that occur in "bunches" , leading to an in- 
crease of the variance by a factor 1 + n/v. 

By definition, a black body has scattering matrix 
5 = 0, because all incident radiation is absorbed. If 
the absorption is not strong enough, some radiation will 
be transmitted or reflected and S will differ from zero. 
Such a "grey body" can still be in thermal equilibrium, 
but the statistics of the photons which its emits will differ 
from the negative-binomial distribution (|7.1|) . A general 
expression for the photon statistics of grey-body radia- 
tion in terms of the scattering matrix was derived in Ref. 
( Bee98j ). The expression is simplest in terms of the gen- 
erating function 



n=0 



from which P(n) can be reconstructed via 



P ( n ) = lim ±^e F( V. 



(7.3) 



(7.4) 



F(0 



tSuj 
"27 



InDet [l - (1 - SS*)Zf] . 



(7.5) 



If the grey body is a chaotic resonator, RMT can be 
used to determine the sample-to-sample statistics of S 
and thus of the photocount distribution. What is needed 
is the distribution of the socalled "scattering strengths" 
<7i, era , . . . o~n > which are the eigenvalues of the matrix 
product SS^ . All a n 's are equal to zero for a black body 
and equal to unity in the absence of absorption. The 
distribution function P({a n }) is known exactly for weak 
absorption (Laguerr e ortho gonal ensemble) and for a few 
small values of N In the large- N limit, the 

eigenvalue d ensity p(a) = (En^( cr — * s known in 
closed- form (Bee99), which makes it possible to compute 
the ensemble average of arbitrary moments of P(n). 

The first two moments are given by 

1 N 1 N 

n = vf-^^{^-Vn), Varn = n + i// 2 — ^(l-cr„) 2 . 



n=l 



n=l 



For comparison with black-body radiation we parameter- 
ize the variance in terms o f the effective number u e g of 
degrees of freedom (|Man95T ). 



Varn = n(l + n/v e g), 



(7.7) 



with t'eft = v for a black body. Eq. (|7.6p implies a reduced 
number of degrees of freedom for grey-body radiation, 



" ^E„(i-^) 2 " ' 

Note that the reduction occurs only for N > 1. 
The ensemble average for N ^> 1 is 

iW"= (l + 7) 2 (7 2 + 2 7 + 2)- 1 , 



(7.8) 



(7.9) 



with 7 = UTdwcii the product of the absorption rate a 
and the mean dwell time Tdwcii = 2ir/N5 of a photon in 
the cavity in the absence of absorption. (The cavity has 
a mean spacing 6 of eigenfrequencies.) As shown in Fig. 
ITTl (red solid curve) , weak absorption reduces i/ e Q by up 
to a factor of two relative to the black-body value. 

So far we have discussed thermal emission from ab- 
sorbing systems. The general formula (|7.5p can also be 
applied to amplified spontaneous emission, produced by 
a population inversion of the atomic levels in the cavity. 
The factor / now describes the degree of population in- 
version of a two- level system, with / = — 1 for complete 
inversion (empty lower level, filled upper level). The scat- 
tering strengths a n for an amplifying system are > 1, 
and in fact one can show that a n i— ► l/cr n upon changing 
a i ► —a (absorption rate i— > amplification rate). As a 
consequence, Eq. (|7.9[) can also be applied to an amplify- 
ing cavity, if we change 7 1— *■ —7. The result (blue dashed 
curve in Fig. IT7|) is that the ratio v c g / v decreases with 
increasing 7 = |<j|rcjweii — vanishing at 7 = 1. This is 
the laser threshold, which we discuss next. 
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FIG. 17 Effective number of degrees of freedom as a function 
of normalized absorption or amplification rate in a chaotic 
cavity (inset). The black-body limit for absorbing systems 
(red, solid line) and the laser threshold for amplifying systems 
(blue , dashe d line) are indicated by arrows. Adapted from 
Ref. l)Bee9Sft . 



B. RMT of a chaotic laser cavity 

Causality requires that the scattering matrix S(u>) has 
all its poles f2 m — iT m /2 in the lower half of the complex 
frequency plane. Amplification with rate a > adds a 
term ia/2 to the poles, shifting them upwards towards 
the real axis. The laser threshold is reached when the de- 
cay rate To of the pole closest to the real axis (the "lasing 
mode") equals the amplification rate a. For a > Tq the 
loss of radiation from the cavity is less than the gain 
due to stimulated emission, so the cavity will emit radia- 
tion in a narrow frequency band width around the lasing 
mode. If the cavity has chaotic dynamics, the ensem- 
ble averaged properties of the laser can be described by 
RMT. 8 

For this purpose, we include amplification in the Wei- 
denmiiller formula (|3.3p . which takes the form 



S(u>) = 1 - 2niW j (uj - ia/2 - Uy x W. 



(7.10) 



The poles of the scattering matrix are the complex eigen- 
values of the M x M matrix 



H = H-iirWW^ = J7diag(f2i-iri 



M- 



-iT M )U-\ 
(7.11) 



The statistical properties of a chaotic laser cavity are closely 
related to those of socalled rando m lasers (see Ref. llCao05f) for 
a review of experiments and Ref. |Tur08) for a recent theory). 
The confinement in a random laser is not produced by a cavity, 
but presumably by disorder and the resulting wave lo calization. 
(Alternative mechanisms are reviewed in Ref. jZai09t) .1 



constructed from the Hamiltonian H of the closed cav- 
ity and the M x N coupling matrix W to the outside. 
Because TC is not Hermitian, the matrix U which diag- 
onalizes TL is not unitary. In the RMT description one 
takes a Gaussian ensemble for H and a non-random W, 
and seeks the distribution of eigenvalues and eigenvectors 
of 7i. This is a difficult problem, but most of th e result s 
needed for the application to a laser are known (|Fvo03r ). 

The first question to ask, is at which frequencies the 
laser will radiate. There can be more than a single lasing 
mode, when more than a single pole has crossed the real 
axis. The st atistics of the laser fre quencies has been stud- 
ied in Refs. (|Mis9d[Hac05l : IZai06l ). Only a subset Abasing 
of the N a modes with r m < a becomes a laser mode, be- 
cause of mode competition: If two modes have an appre- 
ciable spatial overlap, the mode which starts lasing first 
will deplete the population inversion before the second 
mode has a chance to be amplified. For weak coupling of 
the modes to the outside, when the wave functions have 
the Porter-Thomas distribution, the average number of 
lasing modes scales as -/Vi as i ng oc Na /3 (lMis98h . 

Once we know the frequency of a lasing mode, we 
would like to know its width. The radiation from a 
laser is characterized by a very narrow line width, lim- 
ited by the vacuum fluctuations of the electromagnetic 
field. The quant um- lim ited linewidth, or Schawlow- 
Townes linewidth (ISch58l) . 



5u = hKTi/I, 



(7.12) 



is proportional to the square of the decay rate To of the 
lasing cavity mode and inversely proportional to the out- 
put power / (in units of photons /s) . This is a lower bound 
for the linewidth when Tq is much less than the linewidth 
of the atomic transition and when the lower level of the 
transition is unoccupied (complete population inversion) . 
While Schawlow and Townes had K = 1, a ppropr i ate for 
a nearly closed cavity, it was later realized (|Pet79l : ISie89h 
that an open cavity has an enhancement factor K > 1 
called the "Petermann factor" . 

The RMT of the Petermann factor was developed in 
Refs. (|Pat0fi iFraOOh . The factor K is related to the 
nonunitary matrix U of right eigenvectors of TL 1 by 



K = (tfUUiU- 1 !!-^) 



on j 



(7.13) 



where the index labels the lasing mode. (In the presence 
of time reversal symmetry, one may choose U^ 1 = U T , 
hence K = [(UU^)oo] 2 ■) If the cavity is weakly coupled 
to the outside, then the matrix U is unitary and K = 1, 
but more generally K > 1. The probability distribution 
P (K I To) of the Petermann factor for a given value of 
the decay rate Tq is very broad and asymmetric, with an 
algebraically decaying tail towards large K . For example, 
in the case A^ = 1 of a single-mode opening of the cavity, 
P(AT|r )cx {K-l)- 2 -^' 2 . 



18 



Acknowledgments 

I have received helpful feedback on this article from Y. 
Alhassid, B. Beri, H. Schomerus, J. Schliemann, A. D. 
Stone, and M. Titov. I am also indebted to the authors 
of the publications from which I reproduced figures, for 
their permission. My research is funded by the Dutch 
Science Foundation NWO/FOM and by the European 
Research Council. 



References 

[AgaOO] O. Agam, I. Aleiner, and A. Larkin, Phys. Rev. Lett. 

85 (2000) 3153. 
[Akk88] E. Akkermans, P.E. Wolf, R. Maynard, and G. 

Maret, J. Physique 49 (1988) 77. 
[Akk07] E. Akkermans and G. Montambaux, Mesoscopic 

Physics of Electrons and Photons, Cambridge University 

Press, Cambridge 2007. 
[AlhOO] Y. Alhassid, Rev. Mod. Phys. 72 (2000) 895. 
[Alh02] Y. Alhassid and T. Rupp, Phys. Rev. Lett. 91 (2003) 

056801. 

[Alt85] B.L. Altshuler, JETP Lett. 41 (1985) 648. 
[Alt86] B.L. Altshuler and B.I. Shklovskii, Sov. Phys. JETP 
64 (1986) 127. 

[Alt96] A. Altland and MR. Zirnbauer, Phys. Rev. Lett. 76 

(1996) 3420; Phys. Rev. B 55 (1997) 1142. 
[Alt97] H. Alt, C. Dembowski, H.-D. Graf, R. Hofferbert, 

H. Rehfeld, A. Richter, R. Schuhmann, and T. Weiland, 

Phys. Rev. Lett. 79 (1997) 1026. 
[Bar94] H.U. Baranger and P.A. Mello, Phys. Rev. Lett. 73 

(1994) 142. 

[Bee91] C.W.J. Beenakker, Phys. Rev. B 44 (1991) 1646. 
[Bee92] C.W.J. Beenakker and M. Biittiker, Phys. Rev. B 46 
(1992) 1889. 

[Bee93a] C.W.J. Beenakker, Phys. Rev. Lett. 70 (1993) 1155. 

[Bee93b] C.W.J. Beenakker and B. Rejaei, Phys. Rev. Lett. 
71 (1993) 3689. 

[Bee94] C.W.J. Beenakker, Phys. Rev. B 50 (1994) 15170. 

[Bee96] C.W.J. Beenakker, J.C.J. Paasschens, and PW. 
Brouwer, Phys. Rev. Lett. 76 (1996) 1368. 

[Bee97] C.W.J. Beenakker, Rev. Mod. Phys. 69 (1997) 731. 

[Bee98] C.W.J. Beenakker, Phys. Rev. Lett. 81 (1998) 1829. 

[Bee99] C.W.J. Beenakker, in Diffuse Waves in Complex Me- 
dia, edited by J.-P. Fouque, NATO Science Series C531, 
Kluwer, Dordrecht 1999 [arXi^qu^ph79808066 . 

[BeeOl] C.W.J. Beenakker and P.W. Brouwer, Physica E 9 
(2001) 463. 

[Bee05] C.W.J. Beenakker, Lect. Notes Phys. 667 (2005) 131. 
[Ber77] M.V. Berry, J. Phys. A 10 (1977) 2083. 
[Ber78] CP. Berman and CM. Zaslavsky, Physica A 91 
(1978) 450. 

[Ber85] M. V. Berry, Proc. R. Soc. London A 400 (1985) 229. 
[Bla91] J.M. Blatt and V.F. Weisskopf, Theoretical Nuclear 

Physics, Dover, New York 1991. 
[Boh84] O. Bohigas, M.-J. Giannoni, and C. Schmit, Phys. 

Rev. Lett. 52 (1984) 1. 
[Blu90] R. Bliimel and U. Smilansky, Phys. Rev. Lett. 64 

(1990) 241. 

[Bro81] T.A. Brody, J. Flores, J.B. French, P.A. Mello, A. 
Pandey, and S.S.M. Wong, Rev. Mod. Phys. 53 (1981) 
385. 



[Bro95] P.W. Brouwer, Phys. Rev. B 51 (1995) 16878. 
[Bro97] P.W. Brouwer, K.M. Frahm, and C.W.J. Beenakker, 

Phys. Rev. Lett. 78 (1997) 4737. 
[Bro03] P.W. Brouwer, Phys. Rev. E 68 (2003) 046205. 
[Bru96] N.A. Bruce and J.T. Chalker, J. Phys. A 29 (1996) 

3761. 

[Bur72] R. Burridge and G. Papanicolaou, Comm. Pure Appl. 

Math. 25 (1972) 715. 
[Cao05] H. Cao, J. Phys. A 38 (2005) 10497. 
[Cas95] M. Caselle, Phys. Rev. Lett. 74 (1995) 2776. 
[Cha93] J.T. Chalker and A.M.S. Macedo, Phys. Rev. Lett. 

71 (1993) 3693. 
[Cha94] A.M. Chang, H.U. Baranger, L.N. Pfeiffer, and K.W. 

West, Phys. Rev. Lett. 73 (1994) 2111. 
[Cha95] I.H. Chan, R.M. Clarke, CM. Marcus, K. Campman, 

and A.C Gossard, Phys. Rev. Lett. 74 (1995) 3876. 
[Chi71] B.V. Chirikov, F.M. Izrailev, and D L. Shepelyansky, 

Physica D 33 (1988) 77. 
[Das69] R. Dashen, S.-K. Ma, and H.J. Bernstein, Phys. Rev. 

187 (1969) 345. 
[DelOl] J. von Delft and D.C Ralph, Phys. Rep. 345 (2001) 

61. 

[Den71] R. Denton, B. Muhlschlegel, and D. J. Scalapino, 

Phys. Rev. Lett. 26 (1971) 707. 
[Dor82] O.N. Dorokhov, JETP Lett. 36 (1982) 318. 
[Dor84] O.N. Dorokhov, Solid State Comm. 51 (1984) 381. 
[Dor92] E. Doron and U. Smilansky, Nucl. Phys. A 545 (1992) 

455. 

[Dys62] F.J. Dyson, J. Math. Phys. 3 (1962) 140. 
[Dys63] F.J. Dyson and M.L. Mehta, J. Math. Phys. 4 (1963) 
' 701. 

[Efe83] K.B. Efetov, Adv. Phys. 32 (1983) 53. 
[Fal94] V.I. Falko and K.B. Efetov, Phys. Rev. Lett. 77 (1996) 
912. 

[FraOO] K. Frahm, H. Schomerus, M. Patra, and C.W.J. 

Beenakker, Europhys. Lett. 49 (2000) 48. 
[Fyo03] Y. V. Fyodorov and H. J. Sommers, J. Phys. A 36 

(2003) 3303. 

[God99] SF. Godijn, S. Moller, H. Buhmann, L.W. 
Molenkamp, and S.A. van Langen, Phys. Rev. Lett. 82 
(1999) 2927. 

[Goo07] J.W. Goodman, Speckle Phenomena in Optics: The- 
ory and Applications, Roberts & Company, Englewood, 
Colorado 2007. 

[Gor65] L.P. Gorkov and CM. Eliashberg, Sov. Phys. JETP 
21 (1965) 940. 

[Got08] K. Gottfried and T.-M. Yan, Quantum Mechanics: 

Fundamentals, Springer, New York 2008. 
[Guh98] T. Guhr, A. Miiller-Groeling, and H.A. Wei- 

denmuller, Phys. Rep. 299 (1998) 189. 
[Hac05] G. Hackenbroich, J. Phys. A 38 (2005) 10537. 
[Hal86] W.P. Halperin, Rev. Mod. Phys. 58 (1986) 533. 
[Hen99] M. Henny, S. Oberholzer, C. Strunk, and C. 

Schonenberger, Phys. Rev. B 59 (1999) 2871. 
[Hua63] L.K. Hua, Harmonic Analysis of Functions of Several 

Complex Variables in the Classical Domains, American 

Mathematical Society, Providence 1965. 
[Hui98] A.C Huibers, SR. Patel, CM. Marcus, P.W. 

Brouwer, C.I. Duruoz, and J.S. Harris, Jr., Phys. Rev. 

Lett. 81 (1998) 1917. 
[Imr86] Y. Imry, Europhys. Lett. 1 (1986) 249. 
[Jal92] R.A. Jalabert, AD. Stone, and Y. Alhassid, Phys. 

Rev. Lett. 68 (1992) 3468. 
[Jal93] R.A. Jalabert, J.-L. Pichard, and C.W.J. Beenakker, 



19 



Europhys. Lett. 24 (1993) 1. 
[Jal94] R.A. Jalabert, J.-L. Pichard and C.W.J. Beenakker, 

Europhys. Lett. 27 (1994) 255. 
[Kim05] Y.-H. Kim, U. Kuhl, H.-J. Stockmann, and P.W. 

Brouwer, Phys. Rev. Lett. 94 (2005) 036804. 
[Kud95] A. Kudrolli, V. Kidambi, and S. Sridhar, Phys. Rev. 

Lett. 75 (1995) 822. 
[Kue08] F. Kuemmeth, K.I. Bolotin, S.-F. Shi, and D.C. 

Ralph, Nano Lett. 8 (2008) 4506. 
[Kuh05] U. Kuhl, H.-J. Stockmann, and R. Weaver, J. Phys. 

A 38 (2005) 10433. 
[Kui09] J. Kuipers, C. Petitjean, D. Waltner, and K. Richter, 

larXiv:0907.2660l 
[LamOl] A. Lamacraft and B.D. Simons, Phys. Rev. B 64 

(2001) 014514.] 

[Lan97] S.A. van Langen, P.W. Brouwer, and C.W.J. 

Beenakker, Phys. Rev. E 55 (1997) 1. 
[Lee85] P.A. Lee and A.D. Stone, Phys. Rev. Lett. 55 (1985) 

1622. 

[Lew91] C.H. Lewenkopf and H.A. Weidenmiiller, Ann. Phys. 

(N.Y.) 212 (1991) 53. 
[Lod98] A. Lodder and Yu.V. Nazarov, Phys. Rev. B 58 

(1998) 5783. 

[Mah69] C. Mahaux and H.A. Weidenmiiller, Shell-Model Ap- 
proach to Nuclear Reactions, North-Holland, Amsterdam 
1969. 

[Man95] L. Mandel and E. Wolf, Optical Coherence and 
Quantum Optics, Cambridge University Press, Cambridge 
1995. 

[Meh67] M.L. Mehta, Random Matrices, Academic Press, 

New York 1991. 
[Mel88a] P.A. Mello, P. Pereyra, and N. Kumar, Ann. Phys. 

(N.Y.) 181 (1988) 290. 
[Mel88b] P.A. Mello, E. Akkermans, and B. Shapiro, Phys. 

Rev. Lett. 61 (1988) 459. 
[Mel89] P.A. Mello and J.-L. Pichard, Phys. Rev. B 40 (1989) 

5276. 

[Mel96] J.A. Melsen, P.W. Brouwer, K.M. Frahm, and C.W.J. 

Beenakker, Europhys. Lett. 35 (1996) 7; Physica Scripta 

T69 (1997) 223. 
[Mis96] T.Sh. Misirpashaev and C.W.J. Beenakker, JETP 

Lett. 64 (1996) 319. 
[Mis98] T.Sh. Misirpashaev and C.W.J. Beenakker, Phys. 

Rev. A 57 (1998) 2041. 
[Mut87] K.A. Muttalib, J.-L. Pichard, and A.D. Stone, Phys. 

Rev. Lett. 59 (1987) 2475. 
[Nag92] K.E. Nagaev, Phys. Lett. A 169 (1992) 103. 
[ObeOl] S. Oberholzer, E.V. Sukhorukov, C. Strunk, C. 

Schonenberger, T. Heinzel, K. Ensslin, and M. Holland, 

Phys. Rev. Lett. 86 (2001) 2114. 



[Obe02] S. Oberholzer, E.V. Sukhorukov, and C. 

Schonenberger, Nature 415 (2002) 765. 
[OstOl] P.M. Ostrovsky, M.A. Skvortsov, and M.V. Feigel- 

man, Phys. Rev. Lett. 87 (2001) 027002. 
[PatOO] M. Patra, H. Schomerus, and C.W.J. Beenakker, 

Phys. Rev. A 61 (2000) 23810. 
[Pet 79] K. Petermann, IEEE J. Quantum Electron. 15 (1979) 

566. 

[Pni96] R. Pnini and B. Shapiro, Phys. Rev. E 54 (1996) 
R1032. 

[Por65] C.E. Porter, ed., Statistical Theories of Spectra: Fluc- 
tuations, Academic Press, New York 1965. 

[Pri93] V.N. Prigodin, K.B. Efetov, and S. Iida, Phys. Rev. 
Lett. 71 (1993) 1230. 

[Sch58] A.L. Schawlow and C.H. Townes, Phys. Rev. 112 
(1958) 1940. 

[Sch99] H. Schomerus and C.W.J. Beenakker, Phys. Rev. 
Lett. 82 (1999) 2951. 

[Sie89] A.E. Siegman, Phys. Rev. A 39 (1989) 1253. 

[Sil03] PC. Silvestrov, M.C. Goorden, and C.W.J. 
Beenakker, Phys. Rev. B 67 (2003) 241301(R). 

[Smi60] F.T. Smith, Phys. Rev. 118 (1960) 349. 

[Ste96] A.H. Steinbach, J.M. Martinis, and M.H. Devoret, 
Phys. Rev. Lett. 76 (1996) 3806. 

[Sto91] A.D. Stone, P.A. Mello, K.A. Muttalib, and J.-L. 
Pichard, in: Mesoscopic Phenomena in Solids, ed. by B.L. 
Altshuler, P.A. Lee, and R.A. Webb, North-Holland, Am- 
sterdam 1991. 

[Sto99] H.J. Stockmann, Quantum Chaos: An Introduction, 
Cambridge University Press, Cambridge 1999. 

[Tur08] H. E. Tiireci, L. Ge, S. Rotter, and A. D. Stone, Sci- 
ence 320 (2008) 643 (2008) 

[Tra94] C.A. Tracy and H. Widom, Commun. Math. Phys. 
159 (1994) 151; 177 (1996) 727. 

[VavOl] M.G. Vavilov, P.W. Brouwer, V. Ambegaokar, and 
C.W.J. Beenakker, Phys. Rev. Lett. 86 (2001) 874. 

[Vav03] M.G. Vavilov and A.I. Larkin, Phys. Rev. B 67 
(2003) 115335. 

[Vel08] I.M. Vellekoop and A.P Mosk, Phys. Rev. Lett. 101 
(2008) 120601. 

[Was86] S. Washburn and R.A. Webb, Adv. Phys. 35 (1986) 
375. 

[Wei09] H.A. Weidenmiiller and G.E. Mitchell, Rev. Mod. 

Phys. 81 (2009) 539. 
[Wie95] D.S. Wiersma, M.P. van Albada, and A. Lagendijk, 

Rev. Sci. Instrum. 66 (1995) 5473. 
[Wig55] E.P Wigner, Phys. Rev. 98 (1955) 145. 
[Zai06] O. Zaitsev, Phys. Rev. A 74 (2006) 063803. 
[Zai09] O. Zaitsev and L. Devch. larXiv:0906.3449"1 



